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As an explicit example of HQET, we construct correlation functions containing a heavy-light axial current, 
both in the static approximation and in full QCD. This enables us to investigate the size of finite mass effects. 
As we use the lattice regularisation, we can also check the mass dependence of discretisation errors. 



1. CORRELATION FUNCTIONS 

The axial current of one light and one heavy 
quark flavour is phenomenologically interesting, 
as the decay constant /b of the B-meson is de- 
fined by a matrix element of such a current. This 
matrix element requires a non-perturbative treat- 
ment, i.e. a lattice calculation. Since the bot- 
tom quark is too heavy to be included directly in 
Monte Carlo simulations on today's computers, 
the use of effective theories like HQET is a neces- 
sity. In a recent paper ^ , we tested HQET for a 
special correlation function, and studied the dis- 
cretisation errors connected with the heavy quark 
mass. 

The operators under investigation are the axial 
current 



Aq{x) = Tplix)jQj5'lp2ix) 



(1) 



of two relativistic quark flavours ipi and "02 , and 
the static-light axial current 

Af-\x)^Mxhol5Mx), (2) 
where the static quark field satisfies 

P+V^h = ^h, P+ = i(l+7o), (3) 

and its dynamics is determined by the Eichten- 
HiU action |] 

5h-a4^^h(x)VS^h(x), (4) 

X 

The correlation functions we use are defined in 
a finite space-time volume L x L x L x L with 
Schrodinger functional boundary conditions, i.e. 



• The gauge field is periodic in space. 

• No spatial boundary conditions need to be 
specified for the static quark field. The rel- 
ativistic quark fields are periodic in space 
up to a phase 9, e.g. (with unit vector k) 

■4!i{x + Lk)^e'^i}i{x). (5) 

• Dirichlet boundary conditions in time. 

The space-like link variables on the time bound- 
aries are chosen as Uk{x)\xQ=t),L = 1- Correla- 
tion functions can contain derivatives C,', 
Ch, C(i with respect to the boundary quark fields 
(see for their definition). 

We apply the Symanzik improvement pro- 
gramme to the axial currents, defining im- 
proved operators and Aq*'**'^ as well as to the 



action 



We use the correlation functions 

-6 



fii^o) = -^5](Ai(:r)C2(y)75Ci(z)),(6) 

(CUu)75C2(v) X 



fi = - 



2L^ 
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u,v,y ,z 



xC2(y)75Ci(z)), 



/•stat,!/ \ 
J A (^o) 



Jl 



y ^ 



(A^'^*^^(a;) X 



xCh(y)75Ci(z)), 



E (cU")75c;;(v) 

u,v,y ,z 

xCh(y)75Ci(z))- 



(7) 



(8) 



(9) 



2. TESTING HQET 

Our aim is to compare HQET with full QCD 
in the continuum limit, which requires renormal- 
isation. A renormalised coupling is defined by 
matching to the MS scheme, 



2 



9l + 0{gl) 



(10) 



where 50 is the bare lattice coupling. The quark 
masses are also renormalised by matching to the 
MS scheme. The MS mass of the first 

quark flavour is chosen to be zero, while we intro- 
duce the dimcnsionless parameter 



z — Lm 



(11) 



to parametrise the mass of the second quark 
flavour. Taking the continuum limit means 
a/L ^ {) while keeping z flxed. We form the 
ratios 



Xi{L/a) 



Yi{z,L/a) = 



/1(V2) 



(12) 



(13) 



in which the boundary field renormalisation con- 
stants and a divergence due to the heavy quark 
self energy cancel. These ratios are expanded as 



(14) 



and analogously for Yf. We introduce a renor- 
malisation scheme for the light-light axial current 
by demanding that the renormalised ratio Yca 
satisfies the current algebra relations at 2: = 0. 
The static-light axial current requires a scale de- 
pendent renormalisation, and we define a renor- 
malised ratio ^match(/^) by imposing the match- 
ing condition 

+0{l/z) + 0{{a/Lf). (15) 

The renormalised ratios are expanded in powers 
of the coupling analogous to equation ([ij) . 

At tree level, the continuum limits and 
Y-I^^\z) can be calculated analytically. Y-^'^\z) 
contains terms of the form e~^, which are not 
analytic in 1/z. However, the analytic part can 
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Figure 1. The ratio Y^ (solid curves) and the 
expansion of its analytic part at order 1/z (long 
dashes), 1/z^ (short dashes), and Xjz^ (dotted 
curves). 

be expanded; the results at different orders are 
shown in figure |l|. For z > 2 the expansion pro- 
vides a good approximation, and the non-analytic 
part does not seem to be significant. 
At one-loop level, HQET predicts that 



X. 



(1) _ 

matchV"'2,MS 



("^2.Ms) 



{Br-7oln(z|)}x(°) 
0((a/L)2), (16) 



with 70 = -1/4^2 
can be obtained as the 1/z 



This means that B^*""* 



limit of 



nstat 



iz) 



7oln(z- 



-i3y{rW(z)-X«(L/a)},(17) 
Xj 



where it is understood that the continuum limit 
is taken on the right hand side. B^^^ is shown 
for two different values of the phase 9 in figure ^, 
where fits of the expected form 



Sl'"* + ./i- + ./2-ln(z) + 0(l/z2) 



z z 

are also shown. The extrapolation results are 



(18) 
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Figure 2. in the continuum limit. The 

dotted curves are fits to the filled symbols, of the 
form predicted by HQET. 



^stat 1^ 



= 0.0) = 
= 0.5) = 



-0.f36(3), 
-0.f37(f). 



(19) 
(20) 



This agrees with a previous result where ma- 
trix elements between quark states were used. 
Turning the argument around, this means that 
using the B^^^ value from HQET describes 
the heavy quark limit of QCD, independent from 
the chosen correlation functions. To our knowl- 
edge, this is the first time that this fact has been 
established in an explicit example. 

3. DISCRETISATION ERRORS 

A further interesting point to study is the quark 
mass dependence of discretisation errors. We dis- 
cover that for not too large masses, the discretisa- 
tion errors both at tree-level and at one-loop or- 
der are roughly proportional to {am^ ms)^- Ho"^" 
ever, for am^yig > 1/4, the behaviour of discreti- 
sation errors changes completely. This is illus- 
trated in figure ||, where B^^^ is shown both in 
the continuum limit and at finite lattice spacing 
L/a = 32. A closer look shows that this be- 
haviour originates from a breakdown of 0(a) im- 
provement when the quark mass in lattice units 
becomes too large, as it was previously observed 
in a perturbative analysis of the Schrodinger func- 
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Figure 3. B^''**(z) in the continuum limit 
(squares) and for L/a — 32 (triangles). 

tional coupling | [lO| . This means that extrapola- 
tions in the quark mass should be carried out in 
the continuum limit rather than at finite lattice 
spacing. 
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